An Integrable Version of Burgers Equation in "Magnetohydrodynamics" 
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It is pointed out that for the case of (compressible) magnetohydrodynamics (MHD) with the 
fields Vy{y,t) and Bx{y,t) one can have equations of the Burgers type which are integrable. We 
discuss the solutions. It turns out that the propagation of the non-linear effects is governed by the 
initial velocity (as in Burgers case) as well as by the initial Alfven velocity. Many results previously 
obtained for the Burgers equation can be transferred to the MHD case. We also discuss equipartition 
Vy = ±Bx- It is shown that an initial localized small scale magnetic field will end up in fields moving 
to the left and the right, thus transporting energy from smaller to larger distances. 



The Burgers equation [1] has been much studied. It 
can be applied to a variety of phenomena, see e.g. [2], [3], 
[4]. Although this equation satisfies a number of prop- 
erties which are similar to hydrodynamics, it is known 
to be integrable. Hence the Burgers equation does not 
have the properties characterizing chaotic dynamical sys- 
tems. However, to some extent such properties may be 
simulated by random boundary conditions [5], [6]. Also, 
the long time behavior of decaying solutions of the Burg- 
ers equation with an initial velocity which is homoge- 
neous and Gaussian has been studied and many inter- 
esting properties of the spectrum have been found [7]. 
Also, a general statistical theory of the Stochastic Burg- 
ers equation in the inviscid limit has been developed [8] . 
For a recent review of "Burgulence" we refer to the paper 
by Frisch and Bee [9]. These results make it clear that a 
number of properties of the Burgers equation are highly 
non-trivial. With this in mind we present an integrable 
generalization of Burgers equation. The new equations 
are related to magnetohydrodynamics (MHD) in a way 
which is analogous to the relation between the Burgers 
equation and the Navier-Stokes equations. Essentially all 
properties found for the Burgers equation can be applied 
to the new equations, but they also contain some new 
features. 

MHD in 1+1 dimension has been considered by 
Thomas [10] long time ago. The fields are the velocity 
Vx{x,t) and the magnetic field Bx{x,t), which satisfy 



dtVx + Vx dxVx = Bx dxBx + v dlv^, 
dtBx + Vx dxBx = Bx dxVx + (1/cr) dxBx 



(1) 



The first of these equations is similar to the Burgers equa- 
tion. Both equations are modeled from the incompress- 
ible MHD equations. However, since incompressibility 
leads to triviality in 1-|-1 dimensions, the equation div 
V = is not enforced, and the total (including the mag- 
netic) pressure is taken to vanish. Also, in these equa- 
tions it is implicit that variation of the density is disre- 
garded. Similarly, the equation div B = is not satisfied. 



It has been shown by Passot [11] that the above equa- 
tions are not integrable. For a review of some of the 
consequences of these equations, we refer to ref. [12]. 

In the following we consider a different form of 1-1-1 
dimensional MHD, where there is really two dimensions, 
X and y, but where the fields are restricted to depend 
only on y. We restrict ourselves to the fields Bx{y, t) and 
Vy{y,t). The equation of motion^ 



+ (vV)v = (V X B) X B - 



then becomes 



dtVy 



VydyVy ^BxOyBx + V d^Vy. 



Similarly, the equation 

9tB V X (v X B) + (1/cr) V^B, 
where a is the conductivity of the fluid, becomes 

dtBx = -dy{vyBx) + (1/(7) dlBx. 



(2) 



(3) 



(4) 



(5) 



The coupled set of equations (3) and (5) can be inter- 
preted by saying that (3) is a Burgers equation with a 
magnetically generated pressure, governed by (5). Like 
in the case of Thomas's 1-f 1 dimensional equations, div 
V = is not satisfied, and we have also disregarded a 
possible variation of the density. Notice however that 
div B = is satisfied in our approach and that the mag- 
netic pressure B^/2 is kept. 

The main difference between our equations and those 
discussed by Thomas is that his equation (1) includes a 
term BxdxVx which stretches the magnetic field lines, and 
which competes with the term VxdxBx which (in higher 
dimensions) breaks or twists the field lines. On the other 
hand, in our case the magnetic field is divergence free, 
and we included the magnetic pressure. 

The reduced model (1) introduced by Thomas is in 
correspondence to special subclasses of of solutions to the 
full MHD equations. In our case the equations (3) and (5) 



*email: polesen@nbi.dk 

^In eqs. (3) and (5) the density p and the vacuum permeability fio should occur. Thus, (V x B) x B on the right hand side 
should be multiplied by 1/ jiop. We assume a constant p. Then the following rescalings "remove" p and po: o"Po o" and 
B ^ B/^pJ^P- 
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are of course derived from the real MHD equations, but 
they do not directly correspond to known subclasses of 
solutions of the MHD equations. However, the reduced 
model proposed in this paper may have some physical 
interest, which we shall discuss at the end of this pa- 
per after having given the analytic solution of eqs. (3) 
and (5). It turns out that the model predicts that mag- 
netic energy which is initially localized at small scales, is 
moved to large distances by the non- linear dynamics. 

Conservation of energy can easily be checked from eqs. 
(3) and (5). With 



and 



dtivy - B^) + -dy{vy " B = V 9^ (v^ - B ^) . (12) 

These two equations are of the Burgers type. 

We remind the reader that the solution of Burgers 
equation 



dtu + udyU = udyU 



is 



uiy,t) = -{y-a{y,t)), 



(13) 



(14) 



one has 
dEtotit) 



dt 



/OO 
dy [vl + Bl) , (6) 
-OO 

/OO 
dy 
-OO 

/OO I _ \2 I /"^ _ \2 

^ ae--^+^^(«)da / e-^i^-^^^^^da. (15) 



where 
a{y,t) 



Assuming no "diffusion at infinity" , i.e. assuming that 
Vy and Bx vanish for y — > ±00, then 



dt 



= -2 J dy (^ly {dyVyf + i {dyB, 



? ■ (8) 



Thus energy is conserved in the limit v,\/ a Q. 

The idea is now to compare the equations (3) and (5) 
to the well known solution of Burgers equation found by 
Hopf [13] and Cole [14], where the diffusive terms in these 
equations are included. We can show that if 



1/a, 



(9) 



then the equations arc intcgrablc. We do not know if 
the equations are still integrable if i/ ^ l/cr, and/or if 
variations of the density p are included according to the 
conservation equation 



9tP + dy{pVy) = 0. 



(10) 



Of course, the full set of equations can be studied numer- 
ically. 

In the following we shall consider the case where (9) is 
satisfied. It should be emphasized that this assumption 
is not supposed to represent a realistic estimate of the 
magnetic Prandtl number = crv, which e.g. for liq- 
uid metals is of the ordcr^ 10^^. Our excuse for having a 
Prandtl number equal to one is primarily that this allows 
a non-trivial solution of eqs. (3) and (5). Furthermore, 
one space dimension is anyhow not realistic. 

Adding and subtracting eqs. (3) and (5) we obtain in 
the special case (9) 

dt{Vy + Bx) + ]-dy{Vy + B^f = V dl{Vy + ) , (ll) 



Here u{y, t = Q) 
we have 



-dytp{y). For the case where = 



u{y,t) = u{b{y,t),t = 0), 
where b{y, t) solves the implicit equation 
h{y,t)=y-tu{b{y,t),{)). 



(16) 



(17) 



Here b{y,t) can be interpreted as a Lagrangian coordi- 
nate, with b{y, 0) = y for t = 0. This solution can 
be obtained by the methods of characteristics or from 
the saddle point in eq. (15) for ^ 0. In this case 
b{y,t) d{y,t). In the case where there are more than 
one solution of eq. (17) for b{y,t) one should take the 
solution which maximizes the expression 



2t 



(18) 



with respect to &, as is obvious from the saddle point 
expansion of (15). 

For the MHD case the solution can be found in terms 
of the initial values 



Vy{y,t = 0) = -dytpyoiy), 
Bx{y,t = 0) = -dyipBoiv)- 



(19) 



The Hopf-Cole solution of Burgers' equation then gives 



vyiy,t) = - 



y - 2 ("■v+B{y,t) + ay-B{y,t)) 



iy,t)), (20) 



^It can be mentioned that in numerical simulations of the Earth's dynamo with Prandtl number 
value of Pm ~ 0.1 (not so far from one!), since this is what is numerically tractable. See e.g. [15]. 



10 one actually uses a 
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and 



"2^ {av+B{y,t) - ay-Biy,t)) 



vO+BO vO-BO 



(21) 



Here 



av±B{y,t) = 



POO 

/ ae~ 


4t/t 


' — oo 




" POO 


(m— 

4t/t 


— C30 





(22) 



and similarly 
V'Wbo(2/.^) 



/oo (j/-o)^ 
«o(«)±^so(a))e"^^ 
-OO 



f 



(23) 



The last forms of eqs. (20) and (21) follow from 



f 

J — C 



d_ 

da 



In case where v is taken to nearly vanish, the resulting 

saddle point simplifies the solutions (20) and (21), and 
the a's are replaced by solutions of the equations 



av±B =y + t {ip'yo{ay±B) ± ■0so(«f±B)) • 



(25) 



These sohitions can also be obtained by solving the orig- 
inal equations without diffusion {i/ = 1/a = 0) hy the 
methods of characteristics. These solutions can be writ- 
ten in a form analogous to eq. (16), 



and 



Bx{y,t) = - {B,ib+{y,t),0) + B^{b-{y,t),0) 

+vy{b+{y,t),0) - vyib-{y,t),0)) , 



Vy{y,-t) = -j^{vy{b+{y,t),Q) + Vy{b-{y,t),Q) 



(26) 



(27) 



+Sx(6+(y,t),0)-Bx(6_(y,t),0)), 

where b±{y,t) solve the equations 

b±{y, t)=y-t Vy{b±{y, t),0)Tt B^{b±{y, t),0). (28) 

Like in Burgers case b±{y,0) = y. Also, there is a sim- 
ple Lagrangian interpretation of eq. (28), since the right 
hand side involves the initial velocity subtracted or added 
to the usual Alfven velocity Bx{b±,0), where it should be 
remembered that the constant 1 / ^/jMyp was absorbed in 
Bx- In the saddle point of eq. (22) for ;/ — » we again 
have ay±B{y,t) — > b±{y,t), with b± given by the domi- 
nant saddle point, as discussed below eq. (17). It is of 



course easy to show directly that eqs. (26) and (27) sat- 
isfy the original equations (3) and (5) with v = 1/a = 0. 

From the solutions (26) and (27) one can read off 
a few simple properties. If the initial magnetic field 
Bx{y,t = 0) vanishes no magnetic field is generated at 
other times. This is already quite obvious from the orig- 
inal equations (3) and (5). To see this result from the 
solution by characteristics we notice that if B^ vanishes 
at i = it follows from (28) that 6+ = 6_, and cq. (26) 
then gives Bx{y,t) = 0. The velocity field will then be- 
have as a solution of the "pure" Burgers equation for Vy. 

A less trivial case is when the initial velocity field van- 
ishes, 



vy{y,t = 0) = 0. 



(29) 



Then we obtain both a magnetic and a velocity field as 
a consequence of the dynamics, namely 

Bx{y,t) = ^{Bx{b+{y,t),0) + Bx{b-{y,t),0)) , 

Vy{y,t) = ^ {Bx{b+{y,t),0) - Bx{b-{y,t),0)) . (30) 



with 



b±{y,t)=yTt Bx{b±{y,t),0). 



(31) 



We see that if the initial magnetic field is constant, no 
velocity field is generated. However, in general a varying 
initial magnetic field is able to generate a velocity field. 

In MHD it has often been discussed whether there is 
equipartition, Vy = ±Bx, after a long time. In realistic 
MHD simulations one does not in all cases find equiparti- 
tion. A recent study [16] finds that in nonhelical hydro- 
magnetic turbulence in the inertial range the magnetic 
energy exceeds the kinetic energy by a factor of 2 to 3. 
The helical case has recently been discussed in reference 
[17]. In our case it follows rather trivially that if the ini- 
tial fields satisfy equipartition, then this will be true for 
all times. In general, we see from eqs. (26) and (27) that 
equipartition in the exact sense requires 

Bx{b.{y,t),0) = vy{b_{y,t),0) and 
Bx {b+{y,t),0), Vy (6+ {y,t), 0) arbitrary, 
or 

Bx{b+{y,t),0)=Vy{b+{y,t),0) and 
Bx{b-{y,t),0),vy{b-{y,t),0) arbitrary, (32) 



for Vy = Bx and Vy 



-Bx, respectively. In the first case 
it follows from eq. (28) that 6_ = y. Hence from the 
first line in (32) it follows that Bx{y, 0) — Vy{y, 0), so the 
initial fields are equal. Thus, equipartition requires very 
special initial fields and is not possible for general initial 
conditions. 

The considerations above do not, however, answer the 
question concerning equipartition after some time has 
passed. This requires the study of (31) for b± after the 
passage of a sufficient time. Let us consider the case 
Vy{y,Q) = and let us take the initial field Bx{y,0) to 
be localized in in a domain D in y. Then eq.(31) shows 
that b± receive non-trivial contributions from b± G D, 
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and these contributions have different signs, i.e. the non- 
trivial part of y— space move to the right and to the left, 
so the original domain D splits up into right and left 
moving domains y £ Dr and y £ Dl- After some time 
has passed these domains have no overlap. Using the re- 
sult (30) we then see that Bx{y,t) also moves to the left 
(right) with value B^{b-, 0)/2 {B^{b-, 0)/2). Thus the 
value of the magnetic field has decreased by a factor of 2. 
At the same time it follows from (30) that the velocity 
has increased from zero to ±Bx{b±)/2. Therefore one 
has equipartition. 

If the initial velocity is non- vanishing it is more difficult 
to estimate the result from eqs. (30) and (31). We shall 
therefore ask what is the tendency after a short time. We 
start by considering initial fields which are proportional, 



B^{y,0) = Xvy{y,0), 



(33) 



where A is some parameter. We now want to solve eqs. 
(26), (27), and (28) perturbatively for small t. Assuming 
that y is not too close to zero, eq. (28) can be solved 
approximately, 



b±{y,t)^y-{l±X)tVy{y,0), 



(34) 



where the assumption that y does not vanish is needed 
in order that the second term on the right-hand side is 
small relative to y. The fields in (26) and (27) can then 
be expanded, using 

Vyib±,0) « vyiy, 0) (1 - t(l ± A) dyVy{y, 0)) , (35) 

to obtain the results 

B,{y,t) « Xvy{y,0){l-2t dyVy{y,0)), 
Vyiy,t)^Vyiy,0){l-t{l + X^) dyVy{y,0)). (36) 

We see that the first order changes in the fields do not 
preserve the initial proportionality. If we start with a 
magnetic field which is much smaller than the velocity, 
A <C 1, then the relative change is larger in B^ than in 
Vy. If, on the other hand A ~ 1, i.e. the initial fields are 
comparable, then the corrections found in (36) arc of the 
same order. If A is large, i.e. the initial magnetic field is 
much larger than the initial velocity, then the correction 
to the velocity is larger than the correction to the mag- 
netic field. Thus it seems that the non-linearity in the 
basic equations tend to increase the smallest initial field, 
for small times. 

We mention that the contributions to the energy from 
the 0{t) corrections in (36) is proportional to 



dy Vy{y,Of dyVy{y,0) = ^Vy{y,Of 



(37) 



which vanishes since Vy{y, 0) must approach zero at infin- 
ity in order that the energy is finite. Thus the corrections 
in (36) give no contribution to the total energy. In this 
argument we have disregarded that the perturbation may 
not be valid very close to y = 0. 

We have also investigated the situation numerically, 
taking the initial fields Vy{y,0) = sin?/ and Bx{y,0) = 



Xvy{y,0), with — 7r/2 < y < tt/2. We find that after a 
long time the fields fluctuate considerably. In general, 
there is no equipartition, except for A = 1. The (fluctu- 
ating) ratio 



R = 



Bl 



+ £2 



(38) 



is maximally of order 0.006 for A = 0.2 after a long time 
{t = 14). For A = 0.9 one gets a maximum R value 
around 0.37. Finally, for A = 1.1 one gets a maximal R 
around 0.5. Of course, if A = 1 the numerical calculations 
give equipartition with i? = 0.5 for all y. 

It is clear that the usual Burgers shock waves are 
present in our case too. From eq. (28) one has that 
the derivatives of b±{y,t) become inflnite for 



t = t{b±) 



1 



v'ib±,0)±B'^{b±,Oy 



(39) 



In general there will actually be more shocks than in the 
"pure" Burgers case, since derivatives of the solutions 
(27) and (26) contain as well as dyb- which are 

inflnite at t{b+) as well as t{b-). 

In our case the usual conservation form of Burgers 
equation are generalized to 



dtVy = dy 



1 2 



Moreover, we also have 



dtB^ = dy -VyBx + - dyBa 
a 



(40) 



(41) 



In the last equation one needs of course to replace l/tr 
by v in order to apply the solutions found in this note. 

As already mentioned, the many highly non-trivial 
properties of the Burgers equation are shared by the so- 
lutions for Vy and B^- This follows simply from the fact 
that eqs. (11) and (12) are of the Burgers type, so any 
property previously derived can be applied to the fields 
Vy + Bx and Vy — B^. Thus, except for accidental cancela- 
tions these properties also hold for the fields Vy and B^. 
To give an example, consider the correlation function 

C{y, t) =< Vy{y, t)vy{0, t) + B^iy, t)S^(0, t) >, (42) 

where we take homogeneous random fields Vy{y,t) and 
Bx{y,t) with ensemble averages. C{y,t) can be obtained 
from the sum of the correlation functions 

< ivyiy,t) + By{y,t)){vyiO,t) + By{0,t)) > (43) 



and 



< {vy(y,t) ~ By(y,t)){vy(0,t) - By{0,t)) > . (44) 

Each of these two correlation functions contain fields that 
are solutions of the Burgers equations (11) and (12). 
Hence we can use known results from Burgulence, for 
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example from references [7] and [9], to obtain informa- 
tion on C{y,t). The total energy spectrum is then given 

by 

Etot{k,t) = — C{y,t) e^'^y dy. (45) 

271" J-oo 

The total kinetic and magnetic energy is then given by 

/oo 
Etot{k,t) dk. 
-OO 

(46) 

By subtracting eqs. (43) and (44) we can also obtain the 
correlator 

D{y, t) =< B^{y, t)vy{0, t) + S^(0, t)vy{y, t) > 

/OO 
F{k,t) e-^'^ydk. (47) 
-OO 

The function F is related to the Lorentz force 

/OO 
F{k,t)dk. (48) 
-OO 

It is now possible to repeat for example the arguments in 
[7] to obtain information on C{y,t) and F{y,t) if we as- 
sume initial fields which are homogeneous and Gaussian 
with initial spectra of the form 

Etot{k,t)=a''\k\"eo{k), and F{k,t) = p''\k\'' foik). 

(49) 

Here n is the spectral index and eo{k) is an even and 
non-negative function with eo(0) = 1 assumed to be even 
and decreasing faster than any power of k at infinity. 
The function fo{k) has similar properties. The energy 
spectrum at times different from zero can now be ana- 
lyzed completely as in [7]. For example, for I < n < 2 
the spectrum Etot{k,t) has three scaling regions. The 
first is for very small k's, where the large eddies are con- 
served and the behavior agrees with the original with 
a time-indopendont constant. The second region is a k^ 
region, and the third region is characterized by a be- 
havior associated with the shocks. The switching 
from the first to the second region occurs for a k value 
aroimd t^i/(2(2-n))^ whereas the shift to the last region 
occurs around l/Vi, except for logarithmic corrections. 
For —1 < n < 1 there is no inner region and the spec- 
trum develops in a self-similar fashion. For a much more 
complete description, we refer to the original paper [7]. 
It would be of interest to see if somewhat similar results 
are valid in higher dimensional MHD. It should also be 
emphasized that eqs. (11) and (12) are the independent 
equations, and hence one has the possibilty to study more 
general situations than those discussed in the previous 
literature. For example, the fields Vy + and Vy — 
may be started out with different random initial fields 
and their spectra will then develop in different ways. We 
have already seen an example of this phenomenon in the 
perturbative calculation. 



Concerning the use of results obtained from the study 
of the Burgers equation we mention that in eq. (3) one 
can add a forcing term / on the right hand side. The 
very interesting study of the forced Burgers equation in 
ref. [8] can then be used in eqs. (11) and (12), which 
would have / on the right hand side. The master equa- 
tion for the probability density functions of Vy -\- B^ and 
Vy — Bx, their differences and gradients, can then be de- 
rived as in [8]. Again we refer to the literature [8] for 
more information. 

We now return to the question as to whether our pro- 
posed equations have some physical relevance, disregard- 
ing the obvious restrictions due to the low dimensional 
structure. Our approach has the property that although 
the fields depend only on one dimension y, the magnetic 
field Bxijj) points in a different direction. Thus we can 
construct an initial state where the magnetic field in the 
a;— direction is localized in y and ask how this field prop- 
agate from the equations of motion. We take the initial 
field to be 

Bx{y, 0) = Bo = const, for — L < y < L, 
Bx{y,0) = otherwise, 

Vy{y,0)=0 for ally, (50) 

which is a small scale localized field if L is not too large. 
This field is a rudimentary version of a flux tube, which 
realistically would be a magnetic field locally pointing e.g. 
in the a;— direction with a cross section in the y, z— plane. 
In our case this cross section degenerates to a line seg- 
ment —L<y<L. Ignoring diffusion, the time depen- 
dent solution can be obtained from eqs. (30) and (31), 

Bx{y,t) = ^Bo, Vy{y,t) = 
for - L + tBo <y < L + tBo, 
Bx{y,t) = ^Bo, Vy{y,t) = -^Bo 
for — L - tBo < y < L - tBo, 

Bx{y,t) = Vy{y,t) = outside these intervals. (51) 

Thus, the initial flux "tube" decays into two new tubes, 
which move towards the left and the right. Also, half 
of the initial magnetic energy is converted to kinetic 
energy, and equipartition is actually obtained. We see 
that the original localized configuration turns into a less 
localized configuration. This example can be general- 
ized to more complicated initial localized fields where Bx 
has different values in different nearby y— intervals. In 
such cases again the fields in each interval start to 
split and move out to larger distances to the left and 
the right. In this way a fairly localized initial magnetic 
field will end up as a rather delocalizcd state, moving 
energy from small to large distances. This is similar to 
the phenomenon of an inverse cascade in higher dimen- 
sions. This effect has been observed in two and three 
dimensions in realistic simulations of MHD. Since this 
phenomenon occurs in our reduced model one can say 
that our model given by (3) and (5) has some physical 
relevance. However, in three (but not two) dimensions 
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the inverse cascade is usually linked to helicity, which 
does not exist in our case for obvious reasons. 

There are of course very important differences between 
MHD in one and in higher dimensions. For example in 
the latter case (differential) rotation is possible and one 
can have the dynamo effect. Further, important topolog- 
ical changes of flux tubes are possible. Also, the higher 
dimensional MHD equations exhibit genuine chaotic be- 
havior [12], which can only be simulated to some extent 
in 1+1 dimension by having random initial fields. 

I thank Axel Brandenburg for telling me about the 
work in reference [12]. I also thank the referee for some 
interesting comments. 
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